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Abstract. We provide a fairly complete discussion of the electronic properties of nanochains by mod- 
elling the simplest quantum nanowires within a recently proposed approach which combines the Exact 
Diagonalization in the Fock space with Ab Initio calculations (EDABI method). In particular, the micro- 
scopic parameters of the second-quantized Hamiltonian are determined, and the evolution of the system 
properties is traced in a systematic manner as a function of the interatomic distance (the lattice parameter, 
R). Both the many-particle ground state and the dynamical correlation functions are discussed within a 
single scheme. The principal physical results show: (i) the evolution of the electron momentum distribution 
and its analysis in terms of the Tomonaga-Luttinger scaling, (ii) the appearance of mixed metallic and 
insulating features {partial localization) for the half-filled band case, (Hi) the appearence of a universal 
renormalized dispersion relation for the electron energy, which incorporates both the band-structure and 
the Hubbard-splitting features in the presence of electron interactions, and (iv) the transformation from 
a highly-conducting nanometallic state to the charge-ordered nanoinsulator in the quarter-filled case. The 
analysis is performed using the Wannier functions composed of an adjustable Gaussian Is-like basis set, 
as well as includes a long-range part of the Coulomb interaction. 

PACS. 73.63.-b Electronic transport in nanoscale materials and structures - 31.15.Ar Ab initio calculations 
- 71.10.Hf Lattice fermion models - 71.27.-|-a Strongly correlated electron systems 

1 Introduction posed |S] to combine the two efforts in an exact manner, 

at least for small systems. 

In our method of approach (EDABI), we first rigor- 
Recent developments in computational as well as analyti- o^sly determine the ground-state energy Eg of the system 
cal methods have lead to a successful determination of the of mteractmg particles usmg the occupation-number rep- 
electronic properties of semiconductors and metals start- resentation, which is expressed as a function of the micro- 
ing from LDA [T], LDA+U 0, and related |2| approaches, ^copic parameters. Second, we optimize this energy with 
Even strongly correlated systems, such as V2O3 (which respect to the wave-functions contained in these parame- 



undergoes the Mott transition) and high-temperature su- ^^rs by deriving the self-adjusted wave equation for them. 

H ; perconductors have been treated in this manner ^. How- Physically, the last step amounts to allowing the single- 

^- ' ever, the discussion of the metal-insulator transition of particle wave functions to relax in the correlated state, 

the Mott-Hubbard type is not yet possible in a systematic 1^ practice, we propose the particular class of those func- 

manner, particularly for low-dimensional systems. These ^^ns, which are obtained by minimizing variationally the 

difficulties are caused by the circumstance where the elec- ground-state energy Eg with respect to their size (effec- 

tron-electron interaction is comparable, if not stronger, ^^e Bohr radius). In brief, our method of solution does 

than the single-particle energy. In effect, the procedure '^o* ^^i* itself to an exact diagonalization of the paramer- 

starting from the single-particle picture (band structure) i^ed Hamiltonian, but also involves an adjustment of the 

and subsequently including the interaction via a local po- single-particle wave function to obtain a true ground state 

tential, may not be appropriate. In such situations, one of a correlated quantum many-body system, 

resorts to parametrized models of correlated electrons, The EDABI method has been overviewed in a num- 

where the single-particle and the interaction- induced as- ber of papers [MZIIlli so here we concentrate on its ap- 

pects of the electronic states are treated on equal footing, plication to one-dimensional (ID) nanochains of iV ^ 16 

The single-particle wave-functions are contained in the atoms, close to the metal-insulator crossover transition, 

formal expressions of the model parameters. We have pro- This paper complements our recent study of such sys- 
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terns iniC] with the systematic analysis of both half- and 
quarter-filled band cases, as well as with the analysis of 
its transport properties. Throughout the paper we use the 
adjustable Wannier functions composed of a Gaussian ba- 
sis set (ST0-3G), which are determined exphcitly from 
minimization of the system ground-state energy Eq as a 
function of the interatomic distance R. 

The question of delocalization of atomic states also 
has practical relevance. Namely, in dealing with electronic 
properties of quantum dots one usually assumes the exis- 
tence of the effective-mass states. Whether this assump- 
tion is well founded for particular systems can only be 
determined by finding the critical interatomic distance, 
above which the states are localized. In this context, we 
find such critical distance for model systems composed of 
s-like states. 

The structure of the paper is as follows. In Sectional 
we briefly present the basic idea of the EDABI method. 
Then, in Section|3 we provide the numerical ground-state 
analysis of nanochains with the long-range Coulomb in- 
teraction. Namely, we first discuss the electron momentum 
statistical distribution function for the half-filled case and 
analyse it in terms of Tomonaga-Luttinger scaling (includ- 
ing the logarithmic correction, cf. Section I3.1|l and illus- 
trate the charge-density wave ordering for the quarter fill- 
ing {cf. Section lT^ . as well as the analysis the A^-particle 
wave function localization according to the method devel- 
oped by Resta jH]. In Section^ we calculate the system 
charge and spin gaps, as well as perform the finite-size 
scaling with 1/N ^ on these quantities. We also ana- 
lyze the spectral density and extract from it the renor- 
malized dispersion relation for the interacting electrons in 
a nanochain. To the best of our knowledge, such an ex- 
act renormalized one-electron band structure has not been 
determined before. Finally, the application of the present 
scheme to the higher (ns-like) valence orbital systems 
would make possible a direct comparison with the experi- 
mental results for quantum nanowires made of noble and 
alkaline elements. Also, the correspondence between the 
localization criteria infered from the ground-state proper- 
ties on one side, and from the dynamical correlation func- 
tions on the other, is established. 



2 The combined exact diagonalization 
ab-initio (EDABI) method 

The basic idea of the EDABI method is illustrated 
on the block diagram exhibited in Fig. ^ We start from 
choosing the initial Wannier basis set {wi{r)}, composed 
of atomic-like wave functions (here a trial Gaussian basis 
set) with radii a^. Next, we write down the system Hamil- 
tonian in the second-quantization form and determine the 
ground-state energy Eq together with the corresponding 
state in the Fock space |!f'o} by employing the Lanczos pro- 
cedure, for example. Than, the Wannier basis set {wi{r)} 
is optimized with respect to the atomic radii , contained 
in the atomic wave functions (here represented by Gaus- 
sians) composing Wi{r), until the minimal ground-state 



Single-particle 
Schrodinger eq. 



Single-particle basis 



Field operators 



Diagonalization 
in the Fock space 



"1 



*(r), #t(r) /"single-particle "\ ^ I {„,p"(r)} 



I basis optimization J 



Lff = |*o>-Bg<*o| + .. 



*"°(i-),('S'=»)t(,-) 



Ground-state energy Eg — {^o\ H \^o) 



*ij™(ri,...,rv) 



Renormalized N-particle 
wavefimction 

Fig. 1. Flowchart diagram of the EDABI method. The part 
on the far right provides renormalized Wannier functions, field 
operators, and A/^-particle wave function. The top line is ab- 
sent if we select the trial Wannier functions as composed of 
adjustable Gaussian orbitals (e.g. the STO-3G basis). 



energy Eq is reached [H] for a given lattice parameter R. 
On the right, we list the renormalized quantities (for the 
optimized basis), which can be calculated once the whole 
two-step procedure has provided convergent results. 



3 The correlated electrons in a nanochain 

We consider the system of electrons on N lattice sites 
arranged periodically, each site containing a single va- 
lence orbital and an infinite-mass ion (i.e. we start from 
hydrogenic-like atoms). The Hamiltonian, including all 
the direct Coulomb-interaction terms and neglecting other 
(e.g. exchange terms), can be written down (up to a con- 
stant) in the form 



H = efY^r^, + tY^ 



ala^^i+la + HC 



i<j 



(1) 



where 6n, = n, - 1, ef ^ ea + N'^ E»<j(2/% + K,^) 
(in Ry) is the effective atomic level, Rij is the distance 
between the «-th and j-th atoms, t = {wi\T \wi+i) is 
the nearest-neighbor hopping, U = {wiWi\V \wiWi) and 
Kij = {wiWjl V \ wiWj) are the intra- and inter-site Cou- 
lomb repulsions, respectively. In the present form of the 
Hamiltonian, all the mean-field Coulomb terms are col- 
lected in e^^, whereas the last term represents the correla- 
tion part of the long-range Coulomb interaction. We shall 
test a posteriori whether the tight -binding approximation 
for the hopping term is valid. Also, the effect of the direct 
(Heisenberg) exchange is negligible, since the kinetic ex- 
change term will always be dominant [J] . The microscopic 
parameters are expressed in terms of the Wannier func- 
tions {wi{r)} composed of Gaussian-type orbitals. 

The Hamiltonian Q is diagonalized in the Fock space 
with the help of Lanczos technique. As the microscopic 
parameters e^^, t, U, and Kij are calculated numerically 
in the Gaussian basis, the orbital size of the Is-like state 



A. Rycerz and J. Spalek: Fundamental properties, localization threshold, and Tomonaga-Luttinger behavior 



3 



expressed in this basis is subsequently adjusted to obtain 
the minimal ground-state energy Eq, as a function of 
the interatomic distance R. Earlier, we have shown j2| 
that such a combined exact diagonalization - ab initio 
study of the one dimensional system provides the local- 
ization threshold, the electron-lattice couplings, and the 
dimerization magnitude. Moreover, the utilization of the 
Gaussian-type orbitals leads to a variational procedure 
that converges rapidly with the lattice size N |7lll0| . In 
effect, one can extrapolate the optimal orbital parameters 
for larger N using those obtained for small systems (i.e. 
for TV = 6 10), which speeds up the computation re- 
markably. The purpose of this paper is to discuss basic 
solid-state properties of the nanoscopic systems per se, 
and (in some instances) their infinite correspondants by 
performing the finite-size scaling. 

3.1 Statistical distribution and the 
Tomonaga-Luttinger scaling: The half-filled case 

We now discuss the electron momentum distribution of 
ID chain of iV = 6 16 atoms to address the question of 
whether the system composes either a Luttinger -liquid or 
forms an insulating (Mott-Hubbard) state. We first sum- 
marize, following Voit ^Jj, the properties of ID conduc- 
tors, which include the two principal characteristics: 

(z) A continuous momentum distribution function, show- 
ing a singularity near the Fermi level k ~f kp of the 
form (Solyom, Ref. |1H ) 

Uka np + A\kp - kf sgn{kF - k), (2) 

where is a non-universal (interaction-dependent) ex- 
ponent; in consequence, it leads to the absence of fer- 
mionic quasi-particles (the quasi-particle residue in 

vanishes as Zk ~ \kp — kf when k kp). In other 
words, the Fermi ridge is absent in this case. 
(m) Similar power-law behavior of all the other physical 
properties, particularly of the single-particle density 
of states, A/'(cj) |a; — fif (i.e. a presence of a pseu- 
dogap), that imphes a Drude weight Z? > for < 1. 

In the case of lattice models, such as the (extended) 
Hubbard model, the Luttinger liquid behavior is predicted 
by the renormalization group (RG) mapping onto the To- 
monaga-Luttinger model |11| . Through such mapping, one 
can also expect, with the increasing N, a gradual conver- 
gence of the discrete momentum distribution Uka into the 
continuous power-law form ((21 . This hypothesis was first 
checked numerically for the Hubbard model by Sorella et 
al. [12 . 

Here we present the approach to a finite ID chain 
with a long-range Coulomb interaction, as described by 
the Hamiltonian with a simultaneous evaluation of 
the model parameters by optimizing the single-particle 
wave functions contained in those parameters IQ. The 
values of the parameters in the ground state are shown 
in Table as a function of the interatomic spacing (all 
quantities are in atomic units). Only the value Ki of the 



Table 1 . Microscopic parameters (in Ry) of Hamltonian Q for 
a nanochain calculated in the adjusted STO-3G basis compos- 
ing the Wannier functions. The numerical extrapolation with 
A'^ ^ oo is performed. The values of the inverse orbital size 
ct-mi-a (in the units of the Bohr radius ao) and of the ground- 
state energy Eq/N (for N — 10), are also provided. 



R/ao 




ctt 


t 


U 


Ki 


Eg/N 


1.5 


1.363 


0.100 


-0.831 


2.054 


1.165 


-0.749 


2.0 


1.220 


-0.550 


-0.442 


1.733 


0.911 


-0.930 


2.5 


1.122 


-0.797 


-0.264 


1.531 


0.750 


-0.979 


3.0 


1.062 


-0.902 


-0.171 


1.407 


0.639 


-0.991 


4.0 


1.013 


-0.971 


-0.080 


1.291 


0.493 


-0.992 


5.0 


1.004 


-0.987 


-0.037 


1.258 


0.399 


-0.992 



intersite Coulomb interaction Kij = K^i_j^ for nearest 
neighbors is listed, since more distant interactions scale 
essentially in the same manner as their classical values, 
Kij w 2/ |Ri — Rj|. One should also note that the orbital 
size renormalized by the electron-electron interaction 
is about 30 — 40% smaller in the correlated state than 
the corresponding value in the atomic limit (a^^ — qq). 
Therefore, the tight-binding approximation made for the 
hopping term in Eq. is applicable to a good accuracy 
even for the lattice constant R/uq — 2.0, but not much 
below this value. 

The discrete electron-momentum distribution for the 
half-filling {N^ = N) is depicted in Fig. while in Fig. 

we replot it on a log-log scale together with the fitted 
theoretical curves, as explained below. In order to obtain 
these results, we use the boundary conditions (BC) that 
minimize the ground-state energy for a given N (namely, 
the periodic BC for N = -in + 2 atoms and the antiperi- 
odic BC for N = An, at the half-filling). A systematic 
arrangement of the distribution-function data in Fig. [2^ 
for different N values is striking, in that for smaller R it 
is Fermi-like). 

To extract the Luttinger-liquid exponent 9 accurately 
from the data for finite N, it was necessary to also include 
the higher scaling corrections to (O . They can be obtained 
from the Tomonaga mapping in the form of an expansion 
in the powers of hi(TT/\kp — k\R), namely, 

\n.\np -nka\ = -6* In z + felnln z + c + 0(1/ In z), (3) 

where z = TT/\kp — k\R. This singular form of the expan- 
sion is required by the especially slow approach to the 
RG fixed point (Solyom, Ref. [2]). Obviously, by neglect- 
ing the logarithmic corrections one reaches the asymptotic 
form for k ^ kp. The solid lines in Fig. represent 
the formula Q; the fitted values of the parameters 9, b, 
and c are also listed in Table|21 The quality of the fit is de- 
cesively worse for points far away from the Fermi momen- 
tum, and depends on TV since the Fermi wave vector is A'^- 
dependent, i.e. k^ = kf{l - 2/N), where kf = tt/{2R) 
represents the Fermi wave vector in the N —^ oo limit. 
The exponent 9 is also plotted in Fig. [SJ; as a function of 
the lattice parameter R showing that it crosses the critical 
value 9 — 1 (corresponding to the metal-insulator bound- 
ary in ID) for i?c = 2.60ao (ao = 0.529 A is the Bohr 
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(a) 



N=6 
N=8 
N=10 
N=12 
N=14 



-1.0 -0.5 0.0 0.5 

MOMENTUM, kR/TT 

1.6 



1.0 




1.5 2.0 2.5 

ln(jr/|Ap - 



2 3 4 

INTERATOMIC DISTANCE, R (a.u.) 



Fig. 2. Statistical distribution of electrons and Luttinger- 
liquid scaling for a half-filled ID chain of A*' = 6-^14 atoms with 
long-range Coulomb interactions: (a) momentum distribution 
for electrons in the linear and (fe) log-log scale; (c) Tomonaga- 
Luttinger model exponent 9 vs. lattice parameter R (specified 
in ao) and {inset) the corresponding residual sum of squares. 
Solid lines in Figs, (a) and (6) represent the fitting of Eq. @. 



radius). We also give the residual sum of squares {cj. in- 
set in Fig. 12;), which shows that the quality of the fit 
becomes worst for R ^ Rc where the system approaches 
the localization threshold. 

The results for the half-filled system with the on-site 
Hubbard repulsion alone, and with the atomic energy part 
included explicitly as a function of R [T^ , are qualitatively 
very similar to those displayed in Fig.|21 The critical value 
of the lattice parameter in this case is Rc = 2.16ao, and 
does not differ drastically from the previous one (c/. Ta- 
ble 13 for the corresponding values of all the fitted param- 
eters in the Hubbard-model case). This is because such 
nanoscopic systems always have a finite conductivity in 
the large-density limit, since electrons tunnel through a 
finite-width and finite-height potential barrier. Therefore, 
such half-filled systems, both with and without inclusion 
of the long-range interactions, can be considered as be- 
ing close to the metal-insulator transition, in no apparent 
contradiction with the infinite-chain RG result by Fab- 
rizio and the Hubbard-model solution by Lieb and 
Wu |E1 ) which both provide only the Mott insulating be- 
havior. This discussion is complete only after calculating 
the charge and spin gaps, as well as the electric conduc- 
tivity, which are dealt with in the next two Sections. 

The present analysis supplements the earlier discus- 
sion |Kir7lll()| . in which we have interpreted the distribu- 
tion nicer in Fig. in terms of the modified Fermi dis- 
tribution for an almost localized Fermi liquid. The points 



Table 2. The parameters of the expansion J^J for the half-filled 
chain with long-range Coulomb interactions. The correspond- 
ing standard deviations (t(A) for the quantities X = 9,b and 
c are also specified. 



R/ao 


e 


a{e) 


b 


a{b) 


c 


a{c) 


1.5 


0.138 


0.015 


0.147 


0.024 


-0.567 


0.015 


2.0 


0.387 


0.055 


0.425 


0.089 


-0.346 


0.053 


2.5 


0.893 


0.122 


0.971 


0.196 


0.084 


0.118 


3.0 


1.307 


0.128 


1.315 


0.207 


0.357 


0.125 


4.0 


1.455 


0.186 


1.113 


0.299 


-0.032 


0.180 


5.0 


1.413 


0.133 


0.943 


0.214 


-0.823 


0.129 



Table 3. The fitted parameters of the singular expansion Q 
for a half-filled ID Hubbard chain. 



R/ao 


e a{e) 


b 


a{b) 


c 


a{c) 


1.5 


0.229 


0.030 


0.237 


0.048 


-0.537 


0.029 


2.0 


0.803 


0.100 


0.855 


0.162 


-0.078 


0.097 


2.5 


1.283 


0.109 


1.259 


0.176 


0.217 


0.106 


3.0 


1.420 


0.075 


1.230 


0.121 


0.116 


0.073 


4.0 


1.436 


0.069 


1.033 


0.111 


-0.456 


0.067 


5.0 


1.371 


0.037 


0.873 


0.060 


-1.218 


0.036 



are arranged in an almost flat manner for the interatomic 
distance R = 1.5ao, suggesting that some kind of quasi- 
discontinuity of nto- exists near kp. However, an ambiguity 
arises because of the circumstance that for nanosystems, 
the points very close to the Fermi momentum for N —> oo 
system, i.e. the point fc^ — 7r/(2i?), are simply missing. 
Nonetheless, it is amazing that the momentum distribu- 
tion can be rationalized in such simple terms (as either 
the Fermi or the Tomonaga-Luttinger liquids) , which rep- 
resent concepts borrowed from the N oo limit. One 
should also say that the critical value of Rc ~ 2.60ao 
for the localization of the single-particle states obtained 
here is about 30% lower than the corresponding value 
Rc — 3.93ao obtained when we treat the distribution 
nico- as the modified Fermi distribution. However, a dis- 
creapancy of this order should not be suprising anybody, 
since the localization criteria should be treated as semi- 
quantitive at best. One can hope to clarify the situation by 
extending the present analysis to larger N. Nevertheless, 
if we regard a nanoscopic system of A'^ ^ 10 atoms as a real 
systems, then the ambiguity of the statistical distribution 
is significant. Also, the partial localization of electrons in 
a nanosystem will become apparent when we discuss the 
multiparticle wave-function localization at the end of this 
Section, and the conductivity in the next Section. 

We allow ourselves one specific suggestion (if not spec- 
ulation) at this point. Namely, the results shown in Fig. 
13 (and the discussion in the remaining part of the paper) 
point to the possibility that the system with small R < 2ao 
can be analyzed as a Landau-Fermi liquid (albeit with 
discrete momentum states |16jV whereas the system with 
2ao ^ R ^ 2.5ao is closer to the Tomonaga-Luttinger- 
Solyom limit. For R > 2.5ao the electrons can be regarded 
as effectively localized. Such a division into three physically 
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distinct regimes requires a further discussion, carried out, 
to some extent, below. 

Finally, it should be noted that a well-defined spin- 
spin correlations of the spin-density-wave type develop 
with the increasing R for the half-filled discussed 
elsewhere |7ltT^ . The system ground-state energy as a 
function of R is also provided there. 



3.2 Onset of the charge-density wave state for the 
quarter filling 

The electron quasi-momentum distribution for the quar- 
ter-filled (QF) case {N^ = N/2) is shown in Fig. El The 
available number of data points is too small to fit the 
singular formula lO to a reasonable accuracy. In effect, 
the lines on the plots are a guide for the eye. However, 
the smooth behavior of the Luttinger-liquid type is evi- 
dent for R < 4ao, and changes dramatically for the larger 
values of R. This change reflects the onset of the charge- 
density-wave ordering, as illustrated in Fig. ^ In a QF 
chain of = 16 atoms (c/. Fig. the charge is al- 
most uniformly distributed for R < 3ao, but the charge 




0.0 

MOMENTUM, kR/ir 

Fig. 3. Momentum distribution Uka for electrons on a chain of 
A'' = 8-^ 16 atoms in the quarter-filled case {Ne = N/2). Lines 
are drawn as a guide to the eye only. ABC and PBC denote 
the antiperiodic and periodic boundary conditions. 



density wave (CDW) sets in very rapidly in the range 
R/ao — A ~ 5. The CDW order parameter, defined as 
6'cDW = N^^ Z]m(~l)™ {^niAni+,n) , (where Aui = n — 
(n)) approaches its maximal value ^cdw = 1/4 for R > 
8ao (c/. Fig. ^). Also, the crossover range of R, where 
^CDW evolves from 6'cDW ~ to the perfect-order value 
^CDW = l/^i shrinks systematically with the increasing 
N, suggesting quantum critical behavior in the large- A'^ 
limit. One can argue that the charge-ordered state for 
larger R, corresponding to ID Wigner-crystal state on 
a lattice, is unstable in the N ^ oo limit, in accordance 
with the Mermin- Wagner theorem. Namely, we expect the 
amplitude of the quantum charge fluctuations to diverge 
at zero temperature as log A'' for the system with a long- 
range (~ 1/r) Coulomb coupling. However, the divergence 
is absent for R > 5ao in the exact-diagonalization data 
for A^ = 8-^16 (c/. Fig.^), indicating that the zero-point 
charge fluctuations are supressed by the onset of the CDW 
state. In other words, the larger lattice parameter R, the 
larger the system which can be regarded as charge ordered 
in the quarter-filled case. This notion agrees with the ex- 
pected atomic-limit charge-order in the QF case. The cor- 
respondence between the appearance of such charge order 
and the system conductivity is discussed in the next Sec- 
tion, but first, we supplement our analysis of the system 
ground-state properties with the many-body wavefunc- 
tion localization in the framework proposed by Resta [^. 



3.3 Ground-state wavefunction localization 

The approach (9| to the electron localization in the corre- 
lated state is based on the idea that, although the insu- 
lating or the metallic states of matter are usually charac- 
terized by their excitation spectrum, the qualitative dif- 
ference in dc conductivity must also reflect a qualitative 
difference in the organization of electrons in their ground 
state. Such a concept was first emphasized by Kohn in a 
milestone paper ^7], but a complete treatment, related 
to the Berry-phase theory of polarization, was proposed 
over 30 years later by Resta and Sorella (c/. Ref. JHI)- In 
their approach the complex number 



-(a) N=16 (QF) R=4.0 

R=1.5 — I — R=5.0 B 

R=3.0 -X--- R=8.0 - 

f 

i'S'i 




NEIGHBOR DISTANCE, |i-j| 



0.3 



(b) 






N=16 




N=12 




N=8 



2 3 4 5 6 7 8 
INTERATOMIC DISTANCE, R (a.u.) 



Fig. 4. Charge-density distribution for the quarter-filled 
(Ne — N/2) nanochain: (a) density fluctuation correlation 
function {AriiAuj) vs. distance \i — j\, (fe) charge-density wave 
order parameter for the density-density fluctuation {see main 
text for the definition) vs. interatomic distance R. 



^i{2Tr/L)X 



(4) 



(where AT is a many-body position operator defined in ID 
as A" = J2j ^-iid L — NR is the system length) is used to 
discriminate between a localized A^e-particle ground state 
(where I^at^I — > 1 for large systems) and a delocalized 
one, where z^^ vanishes. Namely, the qualitative measure 
of the electron localization is defined as 



1 

A^: 



2 NR"^ 2 
log|zAr| = -— ^log|zjv| (5) 



(where n = Ng/N is the average electron density), whereas 
the phase of the complex number z^^ is related to the 
macroscopic system polarization. One should note at this 
point that the kind of localization described by Eq. ^ is 
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2 3 4 5 6 7 8 

INTERATOMIC DISTANCE, R (a.u.) 



Fig. 5. Localization parameter for the ground-state multipar- 
ticle wavefunction of the half-filled (a) and the quarter-filled 
(b) nanochains. The values of the localization parameter {see 
main text for the definition) arc specified in the units of R. 

clearly not a feature characterizing the individual single- 
particle orbitals; instead, it is a global property of a mul- 
tiparticle ground-state wavefunction as a whole. That is 
because the operator e*^^'^/^-''''- in Eq. cannot be exp- 
resed as a single-particle operator (like X) but rather as 
a genuine many-body operator. 

The numerical evaluation of the complex number z^r^ 
and the resulting {x'^)c within the Lanczos algorithm is 
straightforward, since the operator g 
diagonal in the position representation. The results for 
a ID system with long-range Coulomb interaction are 
shown in Fig. [S] for both the half- and the quarter-filled 
cases (c/. Fig.|5li and[S|D, respectively). To avoid confusion 
by changing the system length L, we plot the localization 
data in units of the lattice parameter R, namely {x'^)c/ , 
instead of {x'^)c- The values of the localization parameter, 
depicted in Fig. E^l, gradually decay with growing R for 
all available numbers of atoms N. However, the dramatic 
change of the decay character takes place around the value 
of the lattice parameter R « 2.5ao (cf. inset in Fig. [31i), 
where the curves for different N join together and form a 
single one for R > 2.5ao. In contrast, for lower values of R, 
the curves are well separated there by showing the diver- 
gence of the localization parameter with N for R < 2ao. 
The results for the Hubbard model (not shown) again does 
not differ qualitatively from those presented in Fig.Et; the 
coalescence point is located near the value R « 2.2ao. The 



positions of the joining points for both the models (with 
and without inclusion of the long-range Coulomb interac- 
tions) are very close to the corresponding critical values 
Rc obtained form the Tomonaga-Luttinger scaling at the 
begining of this Section. 

The numerical results for the quarter-filled case, shown 
in Fig. inb, are of lesser accuracy than those for the half- 
filling. This is because the complex expectation value z^r^, 
defined by Eq. is itself a sum of many complex num- 
bers with different phases, when calculated in the posti- 
tion representation. If the resulting zn^ is close to zero, 
which is the case for QF for all examined values of R; it 
is, in turn, strongly affected by the computer roundoff er- 
rors, which are also amplified by taking the logarithm of 
I^AfJ^ when calculating {x'^)c from Eq. ((SJ. Nevertheless, 
the evolution of the wavefunction localization parameter 
versus R is qualitatively similar to that for the half-filed 
case: the curves for different system sizes N, depicted in 
Fig. Eb, get very close to each other for R > 4. Sag where 
the charge-density wave state is formed, as discussed in 
the preceding subsection. 

The apparent correspondence between the wavefunc- 
tion localization properties and the nature of electron mo- 
mentum distribution, both discussed in this Section, sug- 
gests that a significant reorganization of the ground state 
takes place when the nanochain is close to the localiza- 
tion threshold. These observations are supplemented in 
the next Section with analysis of the system energy gap, 
the spectral function, and the transport properties. 



4 Spectral and transport properties 
4.1 The charge and spin gaps 

For a further verification, of whether the system is quasi- 
metallic or quasi-insulating in the Luttinger-liquid like 
regime, we first perform an extrapolation with 1/N —f 
of the charge-gap defined (for the half-filling) as 

AEc{N) = E^+\N) + E^-\N) - 2E^{N), (6) 

where Eq" {N) is the ground-state energy of the TV-site 
system containing A^g electrons. The corresponding nu- 
merical results are shown in Fig. |S1 where we use again 
the proper boundary conditions {periodic or antiperiodic, 
depending on N) which minimize the ground-state en- 
ergy. The extrapolation of performed using the 
2-nd and the 3-rd order polynomials in {1/N) provides a 
nonzero value of AEq for wide range of the lattice param- 
eter R. Only for the lowest examined value of i? = 1.5ao 
does the gap AEc reaches zero within the extrapolation 
error; for R > l.bao it is clearly nonvanishing. The gap is 
also significantly smaller than the corresponding Hartree- 
Fock value {dotted line) in the regime R < 4.5ao, sug- 
gesting that some kind of reorganization is present in the 
dielectric properties, e.g. a crossover from the Slater- to 
the Mott-type insulator, as discussed for the parametrized 
models ^Hl- This hypothesis is verified by estimating the 
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Fig. 6. Charge-energy gap obtained through the finite-size 
scaling of the results for the chains of A'^ = 4-4-14 atoms. 
The corresponding Hartree-Fock (HF) value of the magnetic 
gap for an infinite system is also drawn for comparison {dotted 
line). The inset exhibits the original data for different values of 
N, used for the scaling. Note a qualitative difference between 
the scaled (with l/N — > 0) and the actual gap values. 
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Fig. 7. Charge-energy gap obtained through the finite-size 
scaling of the results for the chains of A'^ = 4 16 atoms in the 
quarter-filled case. The 2-nd order polynomial has been fitted 
to perform the extrapolation with 1/A — > 0. Inset provides the 
original data used for the scaling. 
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INTERATOMIC DISTANCE, R (a.u.) R (a.u.) 

Fig. 8. Spin gap for the half-filled (a) and the quarter-filled 
(b) nanochains. The results obtained through the finite-size 
scaling are also shown with the corresponding errorbars. 

point is determined later in this Section, where we calcu- 
late the system Drude weight. 

For the sake of completness, in Fig. |Hlwe provide the 
values of the system spin gap, which is defined as 

AEs{N) = E^'-\N)^E^'=°iN), (7) 

where E^- (N) denotes the lowest eigenenergy of iV~site 
system in the subspace with a given total z-component of 
spin Sz- Through the finite-size scaling with 1/7V — > 0, we 
obtain the spin gap AEs « for both the half- and the 
quarter-filled cases (c/. Figs.lHt and|Hb, respectively), and 
in the whole examined range of the interatomic distance 
R. These results clearly show that the insulating phase of 
the system is of the Mott type in the large- limit. No 
indications of the Slater-type phase, for which AEc = 
As > 0, was found. However, for a nanosystem containing 
~ 10 atoms, one can note that the so-called correlated 
insulator {AEc > > 0)? existing for the small values 
of R (cf. insets to Figs. El and and Fig. IHJ gradually 
transforms into the Mott insulator with increasing R. This 
evolution should be contrasted with the evolution of bulk 
3D systems where the Slater antiferromagnet evolves into 
the Mott insulator {cf. Korbel et ai, Ref. [TH|). 



spin gap below. One should also note a qualitative differ- 
ence between the behavior of AEc for nanoscopic chains 
{cf. the inset) and of their N oo correspondants. There- 
fore, the finite-A^ system contains physically different dy- 
namic properties. 

The situation changes when we consider the quarter- 
filled case, Ne = N/2 {cf. Fig. [71). The paraboHc extrapo- 
lation with 1 /N — > now provides the value of the charge- 
gap AEc ~ (within the error bars) for the lattice pa- 
rameter R < 2ao. In the range of R/ao = 2.5 4.5 the 
gap develops (it is significantly greater than the corre- 
sponding error bars), but a random dispersion of the data 
points suggests an inaccuracy of the extrapolation due to 
the nonanalytic behavior of AEc when the system ap- 
proaches the localization threshold. For R > 4.5ao the 
gap grows smoothly to a limiting value corresponding to 
that for the insulating charge-density wave state, identi- 
fied earlier. The more precise position of the localization 



4.2 Spectrum of single— particle excitations: 
renormalized bands vs. Hubbard subbands 

The evolution of the single-particle spectral density func- 
tion Ak(w) with increasing R is shown in Fig. Eli for the 
half-filled {N^ — N) nanochain described by the Hamilto- 
nian The spectral function is defined in the standard 
manner, namely 

71 

(8) 

where the upper (lower) sign correspond to to > ^ {uj < fi), 
respectively, the n-th eigenstate of the system con- 

taining N particles, E^ is the corresponding eigenenergy, 

and the matrix element {'I'n I'^wal^o^ ) ^ with cj^^ = aj^^. 
and Cj^^ = Okcr, is calculated within the Lanczos tech- 
nique set up by Dagotto j2()| . For plotting purposes, we 
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have used the analytical representation of the Dirac delta 
function S{x) {l/TT)e/{x'^ + e^) with e = 0.01 Ry; this 
numerical trick leads to peaks with nonzero width, which 
otherwise would be discrete. In the nanometallic range 
(i? < 2.5ao), the quasiparticle peaks are well defined, but 
incoherent tails are always present and grow in strength 
with increasing R. In effect, in the intermediate regime of 
R ^ 3ao the lower and upper Hubbard bands are formed, 
which, in turn, continuously evolve into discrete atomic 
levels located at the positions w = and to = Ca + U , 
when R —>■ oo {cf. Fig. These limiting peak positions 
correspond to the ground (H") and excited (II~) atomic 
levels. Possibly the most interesting feature of this spec- 
trum is its incoherent nature for R ~ Sag, where the band 
and the interaction energies are comparable and where the 
Luttinger liquid exponent crosses the critical value 9—1 
{cf. Section IXT)l . corresponding to the localization thresh- 
old. 

A direct picture of the spectrum evolution is provided 
by the renormalized dispersion relation, which is obtai- 
ned from the spectral function by extracting the major 
quasiparticle peak position for each momentum fc, and 
is plotted in Figs. EJ^-d. We can again easily identify: 
(i) the nanometalic range {cf. Fig. for R — l.Soo), 
for which the charge gap is of the same order as the en- 
ergy discretization due to the geometrical quantization of 
the quasiparticle momenta, (m) the intermediate regime 
{cf. Fig. Et, R = 3ao) where the gap becomes signifi- 
cantly wider, and {Hi) the atomic limit {cf. Fig. |^, R — 
Sap) in which two dispersionless manifolds located at en- 
ergies uj = ea and uj = Sa + U appear. One should also 
note that the corresponding results for the antiferromag- 
netic Hartree-Fock solution {solid lines) closely match 







(a) 


R/aQ = l.o 












i£/ao = 8.0 




-1 — . 








MOMENTUM, tH/jr 



Fig. 9. Spectral functions A\^{bj) for the nanochain of TV = 10 
atoms (a) and the corresponding renormalized dispersion re- 
lations (band energies with a gap) for the lattice parameter 
K = 1.5ao (&), R = 3ao (c), and R = Sao (d). The quasiparti- 
cle energies for antiferromagnetic Hartree-Fock solution {solid 
lines) are also shown for the comparison. 



those obtained within the EDABI method, both for the 
atomic and for the nanometalic limits (particularly for the 
upper Hubbard band in the latter case). The electron- 
correlation effects, in contrast, clearly appear (near the 
Fermi momenta) for the intermediate range of R, where 
the Hartree-Fock charge gap is significantly larger than 
the exact one {cf. also the preceding Section). Another in- 
teresting feature of the renormalized dispersion relation is 
that the datapoints for different N {= 10 and 12) compose 
a single {universal) renormalized band-like dispersion re- 
lation ek, provided that the proper boundary conditions 
which minimize the ground-state energy for a particular 
N are applied. 

The results for the renormalized band energy ek, to- 
gether with those for the distribution function nko- {cf. Fig. 
Et) J characterize in a fudamental manner the salient fea- 
tures of the electronic states for the quantum nanoliquid. 
Also note, that the gap in the single-particle spectrum is 
always present, that is partially caused by the absence of 
the k states near the Fermi point fc^, and partially by the 
electron-electron interaction. The role of the geometrical 
momentum quantization is dominant in the large-density 
limit {R < 2ao, cf. Fig. 13)), whereas the Coulomb re- 
pulsion determines the single-particle gap in the large-i? 
limit {cf. Fig. Eli). What is remarkable at this point, is that 
the Hartree-Fock aproximation works well in both limits 
of weak and strong electron correlations mentioned above, 
but fails in the intermediate range {R ^ Sap, cf. Fig. IHt). 
One also can say that we observe a magnetic (Slater) 
gap contribution to the renormalized band structure of 
nanoscopic system even though we only observe SDW- 
like correlations, but no spontaneous symmetry breaking. 
This is another fundamental feature of these systems. The 
above results will now be supplemented with the transport 
properties, which are considered next. 



4.3 Drude weight and optical conductivity 

The real part of the optical conductivity at zero temper- 
ature is determined by the corresponding real part of the 
linear response to the applied electric field JUj, and can 
be written as cr(w) = DS{u!) + CTrog(w), where the regular 
part is 



■S {lu - {En ^ Eo)) , (9) 



whereas the Drude weight (the charge stiffness) D is de- 
fined by 



TT 



D = -^{Wo\T\%)--J2 



(10) 



where the kinetic-energy term T is the same as the sec- 
ond term in Eq. and the diamagnetic current opera- 
tor defined as usual: jp = itJ2ja('^]a'^j+ia- ^ HC). Here 
the states j'Z'n) in Eqs. © and (|10|l are the eigenstates 
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Fig. 10. Optical conductivity for nanochains in the half-filled 
case: (a) normalized Drude weight vs. lattice parameter R 
(specified in ao) and its values obtained through the finite size 
scaling (1/A'^ — » 0); (6) regular part of the conductivity, (Trcg(tt;) 
for A'^ = 12 atoms. 
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Fig. 11. Optical conductivity for the quarter-filled chain: (a) 
normalized Drude weight vs. lattice parameter R and its values 
extrapolated to 1/N (the Aitken method has been used 
to estimate the errors); (b) regular part of the conductivity, 
aie^{u!) for A'' = 16 atoms. 



of the Hamiltonian Q corresponding to the eigenener- 
gies En, again with boundary conditions which minimize 
the ground-state energy for a given system size TV. Ma- 
trix elements {'l^n\jp I'fb) are calculated within the Lanc- 
zos method [301. For plotting purposes, we also again use 
the analytic representation of Dirac delta function S{x) 
(l/7r)e/(a;2 + e^), with e = 0.01 Ry. 

For a finite system of N atoms, D is always nonzero 
due to a nonzero tunneling probability through a poten- 
tial barrier of finite width. Because of that, the finite-size 
scaling with 1/A^ — ^ has to be performed on D. Here we 
use, after Gora et al. 21 , the following parabolic extrap- 
olation 



In ID 



N\ 



a + b{l/N) + c{l/Ny 



(11) 



where D'^ denotes the normalized Drude weight D* — 
~{N/tt)D/ {^o\T\^o) for the system of N sites, which 
provides the value in the range ^ D* ^ 1, and thus 
can be regarded as an alternative order parameter for the 
transition to the localized (atomic) state. 

The results for the ID system of = 6 ^ 14 atoms 
in the half-filled chain case are shown in Fig. llOb . The 
values of Dy used for the scaling l|ll|l are listed in Table 
^ together with the resulting and its relative error 
(we stop at the lattice parameter R for which 0, 



within the range of errorbars). We also provide, for com- 
parative purposes, the data for the system with the on-site 
Coulomb interaction only (i.e. for the Hubbard model) in 
Table In both cases, i.e. that with the long-range in- 
teractions (c/. Table and that with the on-site inter- 
action only (c/. Table [5l), the extrapolated Drude weight 
becomes significantly greater than zero (of 2a value) 
only for small lattice parameter, i.e. for R ^ 2.6ao and 
R ^ 2.1ao, respectively. The limiting values match well 
those for which the Luttinger-liquid exponent crosses the 
critical value 9 = \, corresponding to the localization- 
delocalization boundary. The above results suggest again 
the localization onset in these ID systems at half-filling. 
However, the optical conductivity (Treg(w), drawn in Fig. 
llOb . shows the isolated Hubbard peak uj k, U and no 
intraband transitions present in the delocalized state. Be- 
cause of this fact, and also because of the nonzero value 
of the charge-gap for any R (cf. above), one should re- 
gard both the half-filled systems studied here as the Mott 
insulators in the large N limit. Nevertheless, the finite- 
N results show that the conductivity of the nanoscopic 
chain diminishes by two orders of magnitude when the 
corresponding increase of the lattice parameter R is in 
the range of 40 50%. So, one can consider such system 
as undergoing a transformation either from a nanoliquid 
to the localized spin system at the half-filling or from an 
intrinsically partially localized state, since the intraband 
transitions are absent even in the small- i? range. It would 
be very interesting to confirm these results experimentally. 

The situation again becomes completely different at 
the quarter-filling (QF). Namely, the normalized Drude 
weight of the QF systems of A^ = 8-^16 atoms, depicted in 
Fig.lllb. shows a highly-conducting behavior (D*«l) for 
R < 3.5ao, and gradually transforms to zero in the range 
R/ao = 4^5. Also, the regular part of the conductivity 
ai-cgluj) (cf. Fig. Illb ) comprises intraband transitions in 
the quasi-metallic range, particularly near the onset of the 
charge-density wave state (for N — 16). Such a behavior 
provides the model case for the transformation from a 
nanometal in the small-i? range to the charge-ordered 
system {cf. Fig. 0J for larger R. 



5 A brief overview: Novel features 

We have provided a fairly complete description of the 
electronic states in a finite ID chain within the frame- 
work of the EDABI method, which combines the exact 
diagonalization of the many-fermion Hamiltonian in the 
Fock space with a subsequent ab initio readjustment of 
the single-particle (Wannier) functions. Both the ground- 
state and dynamical properties have been obtained as 
a function of the variable lattice parameter (interatomic 
spacing) R. Our approach thus extends the current theo- 
retical treatments of band and strongly correlated systems 
within the parametrized (second-quantized) models, by 
determining those parameters and, in turn, analyzing the 
correlated state explicitly as a function of R. The single- 
particle wave function is allowed to readjust in the corre- 
lated state within the EDABI method, thus unifiying the 
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Table 4. Normalized Drude weight the extrapolated value DJ^, and its relative error for ID half-filled system with 

long-range Coulomb interaction. The critical value Rc ~ 2.6ao is determined from the condition — 2cr(_D^). 



R/ao 


DU 


Dh 




Ds 


^6 




<7(i>i,)/Di, 


1.5 


0.9225 


0.9420 


0.9563 


0.9727 


0.9822 


0.8008 


0.019 


2.0 


0.6245 


0.7105 


0.7875 


0.8660 


0.9222 


0.2567 


0.129 


2.5 


0.1839 


0.2812 


0.4109 


0.5813 


0.7526 


0.0087 


0.420 


2.6 


0.1269 


0.2096 


0.3315 


0.5065 


0.7009 


0.0033 


0.489 


2.7 


0.0840 


0.1508 


0.2595 


0.4312 


0.6441 


0.0012 


0.566 


2.8 


0.0536 


0.1049 


0.1972 


0.3586 


0.5836 


0.0004 


0.641 


2.9 


0.0315 


0.0706 


0.1456 


0.2914 


0.5208 


0.0001 


0.920 


3.0 


0.0196 


0.0461 


0.1047 


0.2314 


0.4575 


0.0000 





Table 5. Normalized Drude weight D'^, the extrapolated value -DJo, and its relative error for ID half-filled system described 
by the Hubbard model. The critical value Rc ~ 2.1ao is determined as for the long-range interaction case (c/. Table|lJ. 
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DU 














1.5 


0.6173 


0.6742 


0.7342 


0.7973 


0.8640 


0.3294 


0.070 


2.0 


0.1529 


0.2259 


0.3276 


0.4631 


0.6331 


0.0092 


0.364 


2.1 


0.0991 


0.1598 


0.2527 


0.3884 


0.5731 


0.0032 


0.449 


2.2 


0.0619 


0.1094 


0.1896 


0.3190 


0.5123 


0.0010 


0.540 


2.3 


0.0372 


0.0724 


0.1382 


0.2561 


0.4514 


0.0003 


0.639 


2.4 


0.0218 


0.0467 


0.0985 


0.2018 


0.3927 


0.0001 


0.739 


2.5 


0.0124 


0.0294 


0.0687 


0.1561 


0.3372 


0.0000 



second and the first quantization aspects of the many- 
particle states into a single, fully microscopic scheme. 

We start by analyzing the situation with one valence 
electron per atom (the half-filled case), and we include 
the long-range Coulomb interaction. The Luttinger-liquid 
type of electron momentum distribution suggests a crossover 
transition from the quasi-metaUic to the insulating (spin- 
ordered) state with increasing R (the same is true about 
the system without the long-range interaction 13 , but 
the quasi-metaUic behavior is manifested to a much stronger 
degree when the long-range part of the Coulomb interac- 
tions is included). The finite-size scaling with (1/A^ 0), 
performed on the charge-energy gap shows the insulat- 
ing nature of the ground state for the large N limit, in 
agreement with the renormalization-group results for the 
infinite system with two Fermi points jl4| . Such an ap- 
parently dychotomic nature (localized vs. itinerant) of the 
nanoscopic systems is confirmed by their transport prop- 
erties. On the one hand, the Drude weight is nonzero in 
the small R limit, and the localization threshold agrees 
with those obtained from the Luttinger-liquid exponent, 
but the regular part of the optical conductivity exhibits 
the insulating behavior. This is the reason we coined the 
term: a partially localized quantum nanoliquid. The most 
fundamental features of the electrons as a quantum nano- 
liquid is provided in Figs.|2t and|nb~d, where the Fermi- 
like distribution, as well as the renormalized band struc- 
ture appear for small interatomic spacing. These features 
evolve with R ^ Rc into atomic-like through the regime 
with split Hubhard-like subbands. 



An illustrative example of the nanonscopic system with 
a clear transformation from nanometal to nanoinsulator 
with the charge-density wave order is provided with the 
quater-filled nanochain (including again the long-range 
Coulomb interaction). For that system, the Drude weight 
is reduced gradually from its maximal value to zero, and 
other properties evolve analogously with increasing lat- 
tice parameter R. The intermediate range of i?, where the 
crossover takes place, also shrinks rapidly with increasing 
N , suggesting the existence of a sharp zero-temperature 
transition in the large N limit. 

The above analysis, for both the half- and the quarter- 
filled cases is very sensitive to the choice of boundary con- 
ditions. This problem is widely studied in the existing lit- 
erature and its relation to the spontaneous magnetic flux 
appearing in the mesoscopic rings of iV = 4n atoms j221 
has been established before |23| . 

The EDABI method implemented here can also be ap- 
plied to discuss the coupling of electrons to the lattice and 
the dimerization (for a discussion, see the third paper in 
Ref. (HI)- Also, one should incorporate the properties of 
the chain with an odd number of electrons (the bound- 
ary conditions will then involve the complex number do- 
main). Furthermore, the dynamics with single holes and 
impurities in the nanochain would make the system more 
realistic as a quantum wire. Finally, the application of the 
present scheme to higher (ns-like) valence orbital systems 
would make possible a direct comparison with experimen- 
tal results on quantum nanowires made of noble and al- 
kaline elements. Nonetheless, we believe that our present 
analysis models the fundamental features of such quantum 
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nanowires, particularly their quantum-nanoliquid aspects 
(evolving Fermi-like distribution, renormalized bands with 
a gap even in the absence of a long-range order). We 
should be able to see a progress along these lines soon. 
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